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Abstract 

We present the canonical and quantum cosmological investigation of a spatially flat, 
four-dimensional Friedmann-Robertson- Walker (FRW) model that is derived from the 
M-theory effective action obtained originally by Billyard, Coley, Lidsey and Nilsson 
(BCLN). The analysis makes use of two sets of canonical variables, the Shanmugadhasan 
gauge invariant canonical variables and the "hybrid" variables which diagonalise the 
Hamiltonian. We find the observables and discuss in detail the phase space of the 
classical theory. In particular, a region of the phase space exists that describes a four- 
dimensional FRW spacetime first contracting from a strong coupling regime and then 
expanding to a weak coupling regime, while the internal space ever contracts. We find 
the quantum solutions of the model and obtain the positive norm Hilbert space of states. 
Finally, the correspondence between wave functions and classical solutions is outlined. 
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1 Introduction 



The search for a theory of quantum gravity constitutes one of the uppermost challenges in 
theoretical high energy physics. The need for quantum gravity finds its roots within Einstein 
general relativity. Powerful general theorems [1] imply that our universe must have started 
from an initially singular state with infinite curvature. In such circumstances, where the 
laws of classical physics break down, it is unclear how any boundary conditions necessary 
for a description of a dynamical system could have been imposed at the initial singularity. 
Quantum corrections could induce a modification of classical general relativity and strongly 
influence the evolution of the very early universe [2]. 

In the last two decades superstring theory [3] has emerged as a successful candidate for 
the theory of quantum gravity. In cosmology, most of the modifications to general relativity 
induced by superstring theory are originated by the inclusion of the dilaton, axion and various 
moduli fields, together with higher curvature terms that are present in the low-energy effective 
actions. Each of these novel ingredients do indeed lead to new cosmological solutions. A 
remarkable example is given by the so-called pre-big bang cosmologies [4] that are derived by 
the low-energy effective string action. Different branches of the solution are related by time 
reflection and internal transformations - 0(d, d) and, in particular, scale factor duality - that 
follow from the T-duality property of the full superstring theory [5] . According to the pre-big 
bang scenario the universe evolves from a weak coupled string vacuum state to a radiation- 
dominated and then matter-dominated FRW geometry through a region of strong coupling 
and large curvature. Although the pre-big bang model has not yet proven able to solve all 
its pitfalls, such as the existence of a singular boundary that separates the pre- and post-big 
bang branches [6] , it provides a good starting point to investigate high-energy cosmology. 

Recently, it has become evident that the five consistent, anomaly free, perturbative for- 
mulations of ten-dimensional superstring theories are connected by a web of duality transfor- 
mations [3] and constitute special points of a large, multi-dimensional moduli space of a more 
fundamental (non-perturbative) theory, designated as M-theory [7]. Quite interestingly, it 
happens that another point of the M-theory moduli space corresponds to eleven-dimensional 
supergravity, the low-energy limit of M-theory [7, 8] . Assuming that M-theory is the ultimate 
theory of quantum gravity, it is natural to begin exploring its cosmological implications. Al- 
though our understanding of M-theory is still incomplete, there are hopes that some of the 
obstacles of dilaton driven inflation in string theory could be overcome within the new theory. 
The underlying idea is to investigate the dynamics at the extreme weak- and strong-coupling 
regimes of superstring theory from a M-theory perspective, where the existence of eleven 
dimensions seems mandatory. 

Several different approaches to M-theory cosmology have already been explored in the 
literature [9],[10],[12]-[17]. In the framework of the Horava-Witten model [11], M-theory and 
cosmology have been combined in the works of Lukas, Ovrut and Waldram [9]. A somewhat 
related line of research is the brane world by Randall and Sundrum [10], where our four- 
dimensional universe emerges as the world volume of a three-brane in a higher-dimensional 
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spacetime. From a different point of view, Damour and Henneaux have investigated chaotic 
models [12] and Lu, Maharana, Mukherji and Pope have discussed classical and quantum M- 
theory models with homogeneous graviton, dilaton and antisymmetric tensor field strengths 
[13]. Different classes of cosmological solutions that reduce to solutions of string dilaton gravity 
have been discussed in Ref. [14]- [16]. In particular, a global analysis of four-dimensional 
cosmologies derived from M-theory and type II A superstring theory has been presented by 
BCLN [17]. Using the theory of dynamical systems to determine the qualitative behaviour of 
the solutions, the authors find that form- fields associated with the NS/NS and R/R sectors 
play a rather crucial role in determining the dynamical behaviour of the solutions: the R/R 
fields cause the universe to collapse but the NS/NS fields, such as the axion, have opposite 
effect. Quite interestingly, for flat FRW models the boundary of the classical physical phase 
space is a set of invariant submanifolds, where either the axion field is trivial or the four-form 
field strength is dynamically unimportant. This interplay leads to important consequences, 
as the orbits in the phase space are dominated by the dynamics associated with one, or the 
other, or both invariant submanifolds in sequence, shadowing trajectories in the invariant 
submanifold [17]. 

In this paper we discuss the main scenario introduced by BCLN from a canonical per- 
spective. This approach allows us to analyse in depth the physical properties of the classical 
solutions and to obtain a consistent quantum description of the model. We consider the 
bosonic sector of eleven-dimensional supergravity which consists of a graviton and an anti- 
symmetric three-form potential. The theory is compactified to four dimensions by assuming a 
geometry of the form M 4 x T 6 x S 1 , where T 6 is a six-dimensional torus and M 4 corresponds 
to a spatially flat FRW spacetime. The effective theory in four dimensions bears a dilaton 0, 
a modulus field (3 identifying the internal space, a pseudo-scalar axion field a and a potential 
term induced by a four-form field. [See Eq. (5) below.] A brief derivation of the previous 
steps is presented in Section 2, where the two invariant submanifolds are identified. In Sec- 
tion 3 we analyse the first invariant submanifold, where the four-form field is negligible and 
the axion field dominates. In particular, in Subsection 3.1 we discuss the parameter space of 
the classical theory and in Subsection 3.2 we find the Hilbert space of states of the quantum 
theory. This programme is performed by using two sets of canonical variables, the so-called 
Shanmugadhasan variables [18] (forming a maximal set of gauge invariant variables) and the 
"hybrid" variables that diagonalise the Hamiltonian. In Section 4, we repeat the analysis of 
Section 3 for the invariant manifold with negligible axion field and dominant four-form field. 
Finally, our conclusions are presented in Section 5. Two appendixes contain technical details 
and complete the paper. 



2 M-theory cosmology 

In this section we closely follow Ref. [17] and derive the four-dimensional minisuperspace 
effective actions that will be used to discuss the dynamics of M-theory cosmology. 
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The bosonic sector of eleven-dimensional supergravity action S^ 11 ) is 

1 



S (11) = J d^Xyf^W - ^..^F 



l 
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(1) 



where Oj, bi, Ci 



. . . 10, F ai ... a4 = 4d[ ai A a2 _ a4 ] is the four-form field strength of the an- 
tisymmetric three-form potential -A ai ... 03 , and (?( n ) denotes the determinant of the eleven- 
dimensional metric . Equation (1) describes the low-energy limit of M-theory [7, 8]. 

The four-dimensional effective action is derived from Eq. (1) by a sequence of two Kaluza- 
Klein compactifications, on a circle S 1 with radius = e* 1() / 3 and on an isotropic six-torus 
with volume Vt& = e 6/3 , respectively The eleven-dimensional metric is 

ds* U ) = 9ab ) dX a dX h = e-^'^g^dx^dx" + e^S^dxP] + e 2 * lo/3 dX 2 , (2) 

where fi,u = ... 3 and i,j = 4 ... 9. The first compactification gives the ten-dimensional 
effective action for the massless type IIA superstring [7, 20] 
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where g$ n = (g^, e 2/3 %), and H mnp and F mnpq are the field strengths of the potentials B np 
and A npg , respectively. Assuming that the only non-trivial components of the form fields 
are those associated with M 4 the second compactification gives the effective four- dimensional 
action [17] 



S = J d'x^g- [e-* 4 (i? (4) (^) + (V$ 4 ) 2 - 6 (V/3) ; 



12 



48 



(4) 

where the four-dimensional dilaton field is $4 = $10 — 6(3. Finally, solving the field equation 
for the four-form F^ vXk and dualizing the three-form H' IuX , Eq. (4) can be cast in the form 
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(5) 



where a is the pseudo-scalar axion field dual to the three-form, H^ vX = e^e^^V K a, and 
p\iv\K _ Q e -6/3 e ^u\K Eq Ua tion (5) is our starting point to investigate classical and quantum 
four- dimensional M-theory cosmology. 

Since we are interested in homogeneous and isotropic four-dimensional cosmologies the 
ansatz for the four-dimensional section of the (string frame) metric (2) is 



ds 2 ^ 



g^dx^dx" = -N 2 (t)dt 2 + e 2a{t) dVL zk , N(t) > 
4 



(6) 



where dQsk is the maximally symmetric three-dimensional unit metric with curvature k = 
0, ±1, respectively. By substituting Eq. (6) in Eq. (5) and requiring that the modulus field /3, 
the dilaton $ 4 , and the axion a depend only on t, the action (5) becomes 



S = fdt - (sa 2 - <j> 2 + 6f3 2 + I<j 2 e 2 ^ +< «) + pt (dke' 2 ^ - X -Q 2 e^ 



-0-6/3 



n 

where we have defined the "shifted dilaton" field 

4> = $4 — 3a , 

and the Lagrange multiplier 

H(t) = Ne^ > . 



(7) 



(8) 



(9) 



The dynamics of the action (7) has been discussed qualitatively in Ref. [17]. Here we restrict 
attention to spatially flat spacetimes (k = 0) and discuss in detail two subcases of Eq. (7) 
that turn out to be completely integrable: 



I The invariant submanifold 



Sj = [ d t\- ( 3a 2 - 2 + 6/? 2 + -a 2 e 2( - 3a+ ^ 
J ji \ 2 



(10) 



which is obtained from Eq. (7) by setting Q = and k = 0. This case describes spa- 
tially flat low-energy M-theory cosmology with negligible four-form F^ uXk . The general 
solution for this model (including spatially curved models which are not discussed here) 
was first discussed by Copeland, Lahiri and Wands in Ref. [15]. 



II The invariant submanifold 



SlI = fdt -(3a 2 -<j> 2 + 6f3 2 ) - -/iQ 2 e 3 ^ 
J ji v ' 2 



-4>-6f3 



(11) 



which is obtained from Eq. (7) by setting a = constant and k — 0. This case describes 
spatially flat low-energy M-theory cosmology with trivial axion field. 



The intersection of the two invariant submanifolds I and II further identifies the invariant 
submanifold 



S InII = fdt - (3d 2 - 2 + 6/3 2 ) 



(12) 



In the following two sections we will discuss the classical and quantum dynamics of the in- 
variant submanifolds I and II, respectively. The degenerate case Smu will be outlined at the 
end of Section 4. 
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3 Invariant submanifold I 



Equation (10) can be cast in the canonical form 

Sj = J dt dp a + (fiptp + f3pp + ap a -H 
where the Hamiltonian is 

H = nH, H=± (2pl - 6pl +pl + 12^e- 2 ^)) 



(13) 



(14) 



As is expected for a time-reparametrization invariant system, the total Hamiltonian TL is 
proportional to the non-dynamical variable \x [21, 22]. The latter enforces the constraint 

H = - 2p2_ 6 p2 +p 2 +12p 2 e -2(3a + 0) =o (15) 

The canonical equations of motion are 



a 



Pa 



6 ' 

p a = 3ple'^ 3a+ ^ 



r 2 ' 



h = ple- 2 ^ 



1 12' 



P/3 = 0, 

where the dots represent differentiation w.r.t. gauge parameter 



a = p a e 
Pa = 0, 



-2(3a+0) 



(16) 



r(t) 



to 



H(t')dt' , 



t>t , 



(17) 



where to is an arbitrary constant. Note that since /x is positive defined r(t) is a monotone 
increasing function. 

Different choices of /i correspond to different choices of time. We have the following 
interesting time parameters: 

i) Cosmic proper time t c . The cosmic proper time is defined by the condition N — 1. It is 
obtained by choosing 

dt r 



dr 



which leads to 



t c (r)= f dre-W, 
where we have chosen the arbitrary constant so that t(—oo) = 0. 



(18) 
(19) 



ii) Gauge proper time t g . The gauge proper time is defined by the condition /i — 1. This 
leads to 

t g (r) = r-T . (20) 

This is the natural time parameter to discuss the classical dynamics and the quantization 
of the model. 
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Hi) The BCLN time variable rj. This is defined by the condition 

' ln = e (3«+0-6/3)/2 ? ^ 



dt c 

or, using Eq. (19), 

dV = e (3a-<^-6/3)/2 ^ ^2) 



The BCLN dimensionless time variable t bcln 

dT Rnr „ dS 



o = 0. (23) 



, — , ) ' 'BCLN 'BCLN 

arj dr) 

The choices Hi) and iv) are useful to compare the notations used in this paper to those of Ref. 
[17]. This is done in Appendix A. 



3.1 Classical solutions 



The off-shell solution of the canonical equations is 1 

a = a + ^ln [cosh (k(t - r ))] - £(r - r ) , 

p a = 3/ttanh [k(t — r )] — 6£ , 

= 0o-^ln[cosh(/t(r-r o ))] + 3^(r-r o ), 

Ptf, = /ttanh [k(t - r )] - 6£ , 

= A) + y§(r-r ), 

P/3 = constant , 

o" = o" H tanh [k(t — r )] 

Pa 

p CT = constant , 
where a , <p , j3 , a and r are constants of integration, 

_2 10f 2 , P/3 



(24) 



^-12^ + fg = 2if, «^0, (25) 



1 Here and throughout the section we assume c^O (p a ^ 0). The case k = corresponds to the degenerate 
case S/n/z an d will be discussed at the end of Section 4. 
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and (we choose k > for simplicity) 



3a + 



In 



\Po 



K 



(26) 



On the (physical) shell £ = is a degenerate trivial configuration of the system because it 
implies k — and pp = 0. 

Since the model is integrable we can find a maximal set of gauge invariant observables. The 
system is invariant under reparametrizations of time so we expect six physical gauge invariant 
observables. (Indeed, the system possesses four canonical degrees of freedom and thus the 
general solution is determined by eight constant of motion: the constraint H, the initial value 
of the gauge parameter r, and six other constant quantities which are the observables of the 
system.) Considering the off-shell solution a possible choice is 



U = a + + 3 ^ Pa ar ccosh ( 

6/C Vl^l y 



V = \ (3P<t> ~ Pa) , 



where 



X = f3 



PP 



Y 



a — 



12k 

Pa~P<f> 

2p* 



arccosh 



\Pt\ 



4>+3a 



K = 



w 



Pa- 



Pa, 



(27) 



{Pa-P^f +ple 



-2(3a+<f>) 



1/2 



(28) 



The quantities (27), being gauge invariant, have vanishing Poisson brackets with the Hamil- 
tonian (14) and have been chosen such that the only nonvanishing Poisson brackets are 



[X, W], 



(29) 



The set of gauge invariant quantities (27) can be completed by the Hamiltonian H and by its 
canonically conjugate ([T,H] P = 1) 



T 



arccosh 



K 



\P<t\ 



(30) 



to obtain a maximal set of canonical variables (also called Shanmugadhasan variables [18]). 
Note that the quantity T has been chosen to have vanishing Poisson brackets with the observ- 
ables (27). Since [T,H]p = 1, T transforms linearly under gauge transformations generated 
by the constraint. This will be useful in the following. 

For sake of completeness, let us write the relation between the gauge invariant observables 
and the constants of motion in Eqs. (24). Using Eqs. (14), (27) and (30), after a bit of algebra 
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we find 



1 , \Pa 

2 K 



-\u, (f> = -W-^ + ?-U, (3 = X, a = Y, 
2 2 K 2 



K 



— (4V 2 - W 2 ) + 2H 
1 2 



1/2 



(31) 



Another useful canonical chart is formed by the hybrid variables that diagonalise the 
constraint (15). Although the hybrid variables are not (all) gauge invariant they allow to fix 
a global gauge and quantize exactly the system. The hybrid variables (a, b, c, a) are defined 
by the canonical transformation 



+ 3a , 



Pa = ^ (Pa 



b = y/3(<i> + a) , 
1 



(32) 



Pb 



2^/3 



(3p^ - p a ) , Pc 



2^/3 



Note that a coincides with the four- dimensional dilaton field $ 4 . Using the hybrid variables 
the constraint (15) reads (we have divided by a factor 12) 



2 2 i 2 i 2 -2a n 
Pa ~ Pb + Pc + P* e = • 



(33) 



The gauge invariant canonical variables are related to the hybrid variables by the canonical 
transformation 



U = b + — arccosh 
y/3 K 



X = — = c — — arccosh 
2V3 k 



K 



\Pa\ 



K 



\Po 



Y = a , 

Pa 



T 



■ arccosh 



K 



\Pa\ 



V = V3 Pb , 
W = 2V3p c , 

Z =Pa, 

H= 1 -(PI-Pl+Pl+Ple~ 2a ) 



(34) 



where k = ^p\ + p 2 a e 2a . The canonical transformation above is generated by the generating 
function of the first kind [19] 



Fi(a,6,c,<T;C/,X,y,T) 



2T 



(V3U - bf - {2VSX - cf - arccosh y 1 + e 2a (F - of 



(35) 
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Let us now discuss the behaviour of the classical solution. The on-shell classical solution is 
determined by six physical parameters. From the equations of motion we see that five of them 
(ao, 0o, A), &o, and r ) give initial conditions for the canonical variables and do not influence 
the qualitative behaviour of the solution. Therefore, the qualitative dynamics of the model is 
determined by a two-dimensional parameter space described by two coordinates, for instance 
k and £. Using k and £ as free parameters, from Eq. (25) it follows that pp is (on-shell) 

Pp = ±2v / 3 v /l2 £ 2 - k 2 . (36) 

The sign of pp determines the dynamical behaviour of the internal six-torus space. In fact, 
from the solution of the equations of motion one obtains the scale factor of the internal space 
(we set r = for simplicity) 

R T6 = e p = eP° ■ e^ r/12 . (37) 

Physically, we want compactification at late (cosmic) times. 2 Compactification of the six- 
torus space is achieved for pp < 0. Indeed, for negative values of pp lim T ^ 00 R T % = 0, i.e., the 
internal space shrinks to zero for large values of the gauge proper time, while decompatifying 
for r — > — oo when the strong coupling region of the theory is approached. Since the relation 
between the cosmic proper time and the gauge proper time is monotone [see Eq. (19)] the 
dynamics in r is identical to the dynamics in t c , and the internal space shrinks to zero for 
large values of the cosmic proper time as well. Note that the gauge proper time is defined 
in the interval r e] — oo, oo[ whereas t c is defined only on the half line. The limiting value 
Pp = corresponds to a constant (stable) internal space with radius Rt@ = e@° and is physically 
acceptable for sufficiently large negative values of f3 . In the following we will restrict attention 
to nonpositive values of pp, the extension to pp > being straightforward (see Fig. 1). 

At fixed k we distinguish three different dynamical behaviours of the four-dimensional 
external space according to the value of £: 

i) i < — k/2 (pp < —2k). In this case the external scale factor ever expands while the 
internal scale factor shrinks from infinity to zero. The Hubble parameter is always 
positive and vanishes asymptotically at large times. In particular, for £ = —k/2 the 
external space starts at r = — oo (t c = 0) with a finite nonzero scale factor and vanishing 
Hubble parameter. For £ < —k/2 the external space starts with a vanishing scale factor 
and infinite Hubble parameter, which is always decreasing, r = — oo (t c = 0) is the 
strong coupling region where the coupling constants of the theory, g = exp(0) and 
<7io = exp($io), become infinite. Conversely, r = oo (t c = oo) is the weak region 
coupling where both g and gio vanish, g and g\o are always decreasing. (See Fig. 2.) 

2 A successful physical model ultimately requires that the moduli fields are stabilized as well. Stabilization 
of the internal space does not occur in the models under consideration, where only a fraction of all the 
degrees of freedom present in Eq. (1) are considered, with exception of the (fine-tuned) case pp = (see 
below). Hopefully, the inclusion of more degrees of freedom will provide a mechanism for stabilization of 
extra-dimensions at late times. 
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-K/3 +K/3 




Figure 1: Parameter space for invariant submanifold I. The physical points are represented by the 
two branches of the hyperbola (36). The different colored regions determine on the two branches of 
the hyperbola the different physical cases described in the text. 



o 




Figure 2: Graphical representation of the dynamics of the invariant submanifold I for £ < —k/2 as 
function of the gauge proper time t g and of the cosmic proper time t c . [Here and in the following 
figures we set «o = = A) = for simplicity] The solid and dashed lines describe the evolution 
of the radius of the external spacetime and of the internal space, respectively. The dotted and 
dash-dotted lines represent the evolution of the coupling constants g and g±o, respectively. 
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Figure 3: Graphical representation of the dynamics of the invariant submanifold I for — k/2 < £ < 
— k/3 as function of the gauge proper time t g and of the cosmic proper time t c . The internal scale 
factor and the coupling constant g and gio are ever decreasing and vanish for large times. The 
external scale factor is large both at small and large times. (See Fig. 2 for the key to the plotted 
lines.) 

ii) £ > k/2 (pp < —2k). In this case both the external scale factor and the internal 
scale factor ever contract. The Hubble parameter is always negative and asymptotically 
vanishing at small times (r = — oo). In particular, for £ = k/2 the external space ends 
at r = oo (t c = 0) with a finite nonzero scale factor and vanishing Hubble parameter. 
For £ > k/2 the external space ends with a vanishing scale factor and infinite Hubble 
parameter, g increases from zero at r = — oo (t c = — oo) to infinity at r = oo (t c — 0). 
Conversely, gio decreases from infinity to zero. 

Hi) —k/2 < £ < —k/2\[3 and k/2\[3 < £ < k/2. In this case the external scale factor first 
contracts then expands, bouncing from infinity to infinity. In particular, for 

a) —k/2 < £ < —k/3 (pp < —2k) the internal scale factor shrinks from infinity to 
zero. The Hubble parameter starts with infinite negative value, becomes positive 
and then decreases to zero after having reached a positive maximum, g and g w 
decrease from infinity to zero [g w to a finite nonzero positive value for the limiting 
value £ = —k/3] (see Fig. 3); 

b) —k/3 < £ < —k/2\[3 the internal scale factor shrinks from infinity to zero [—k/3 < 
£ < -k/2^3 (-2k <pfs < 0)] or is constant [f = -k/2^3 (pp = 0)].The Hubble 
parameter starts with infinite negative value, becomes positive and then decreases 
to zero after having reached a positive maximum, g decreases from infinity to zero. 
gw bounces from infinity to infinity via a positive minimum; 

c) k/2\[3 < £ < k/3 (—2k < pp < 0). The internal scale factor shrinks from infinity 
to zero \k/2\[3 < £ < k/3 (—2k < pp < 0)] or is constant [£ = k/2\[3 (pp = 0)]. 
The Hubble parameter is first negative and small, decreases to a negative minimum 
and then increases to infinity, g increases from zero to infinity, gio bounces from 
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infinity to infinity via a positive minimum. 

d) /c/3 < £ < k/2 (p/3 < —2k). The internal scale factor shrinks from infinity to zero. 
The Hubble parameter is first negative and small, decreases to a negative minimum 
and then increases to infinity, g increases from zero to infinity. g w decreases from 
infinity to zero for /c/3 < £ < k/2 or to a finite nonzero positive value for the 
limiting value £ = /c/3. 

Several comments are in order about the dynamics of the dilaton-axion-moduli solution 
First of all, we do not find a two-branch solution as happens in the pre-big bang 
scenario [4]. Indeed, the action (10) (with a ^ 0) is not invariant under scale factor duality 
and moreover the presence of the axion leads to an effective potential [see Eq. (33)] with 
opposite sign to that required for the pre-/post-big bang branches to exist. (See e.g. Gasperini, 
Maharana and Veneziano in Ref. [6].) Scenarios i) and Hi) are most attractive as far as 
a physical description of our universe is concerned. According to i) the universe begins in 
a strong coupling region with large internal space. Actually, both the coupling constants, 
g and g w , and the radius of the internal six-torus are infinite at t c = 0. However, this 
is not disturbing because for t c — > the spacetime curvature blows up. So we are in the 
strong coupling regime of the theory and the low-energy description of M-theory given by 
eleven- dimensional supergravity breaks down. Hopefully, non-perturbative effects would cure 
the initial singularity. At small times (t c ~ 0) the expansion of the external spacetime is 
characterized H » 1. Then the expansion slows down and at large times the Hubble 
parameter vanishes, approaching a standard FRW behaviour. At small cosmic times the 
scale factor of the external spacetime expands as 

a(t c ) » 42C/k+i)/(6*/k+i) . (38) 

This behaviour coincides with that found by Copeland, Lahiri and Wands in Ref. [15]. (See 
Eq. (4.43) and (4.46) with £//c = ±< a /A.) For £ < -k/2 the exponent in Eq. (38) takes 
values in the interval ]0, l/3[. This implies that we have no inflationary expansion in the 
model. (Here and throughout the paper we follow Ref. [15] and define inflation as accelerated 
growth of the external scale factor, namely d 2 a/dt 2 c > 0.) Inflation happens instead in case 
Hi), where the external spacetime is first contracting and eventually expanding. For this 
region of the parameter space the internal space is ever decreasing. At small times both 
the internal space and the ten-dimensional coupling constant are decompactified. At large 
times the eleventh dimension may be compactified [ )-d)] or decompactified [cases b)- 

c)]. Among the former ) (Fig. 3) is the most relevant from the physical point of view 

since the Hubble parameter vanish asymptotically at large times. 

In the next subsection we will quantize the model both in the gauge invariant and hybrid 
variables and prove the equivalence between the quantum theories in the two representations. 
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3.2 Quantization 



Since the model is completely integrable the quantization can be performed exactly. Further- 
more, since we have reduced the system to a maximal set of gauge invariant variables, the 
quantization in the Shanmugadhasan representation is straightforward and the reduced and 
the Dirac approach are equivalent [23]. We will discuss the reduced method for the maxi- 
mal gauge invariant representation and the Dirac method for the hybrid representation for 
simplicity. 3 

Let us consider first the gauge invariant canonical representation. Since T transforms 
linearly under gauge transformations, a natural gauge fixing is 

T = t-t . (39) 

By imposing Eq. (39) the Lagrange multiplier is 

H = l, (40) 

so we are working in the gauge proper time. The effective action is 

S cS = Jdt [UV + XW + YZ- H cS ] , (41) 

where 

H cS = -H = 0. (42) 

Let us quantize the model. In order to do this we must choose a representation for the 
quantum observables. Taking into account Eqs. (27) we choose the representation in which 
V, W and Z are differential operators with eigenvalues v, w and z, respectively. The Hilbert 
space measure is 

d[w] = dudxdy , (43) 



and V, X and Z read 



v=-i-?-, ^ =-*!-> ^ =-*!-■ ^ 

ou ox ay 

The Schrodinger equation reads 

d 

H eS y(u,x,y;t) = i—y(u,x,y;t). (45) 

Since the effective Hamiltonian of the system is vanishing identically the wave functions do not 
depend on t, i.e., ^f(u, x, y; t) = ^(u, x, y). (This equation is nothing else that the Wheeler-De 

3 See also Ref. [13] where the quantization is achieved by mapping the interacting classical equations into 
sets of free equations and then introducing the corrisponding intertwining operators. 
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Witt (WDW) equation of the system.) An orthonormal basis in the Hilbert space is given by 
the set of eigenf unctions of the observables (44), namely 



^^^y) = j^r^ l{vu+wx+zy) - (46) 

Let us briefly discuss the physical interpretation of Eq. (46) . From the definition of W in Eqs. 
(27) we find w = pp. Therefore, the eigenstates (46) with w > (w < 0) corresponds to a 
growing (decreasing) scale factor of the internal space. Analogously, z > (z < 0) corresponds 
to growing (decreasing) axion field. From Eqs. (31) we can define the operator corresponding 

tO K 2 

I ( d 2 d 2 \ 

ft2 = i2 (" 4 ^ + a?) ' (47) 

The eigenvalues of the operators k 2 and of the operators in Eq. (44) are related by the relation 
12k 2 = 4v 2 -w 2 . [See also Eq. (36).] Using the previous results it is straightforward to identify 
the wave functions that correspond to the classical cases 

Now let us turn to the hybrid canonical chart. From the constraint (33) it is natural to 
choose 

d[w] = dadbdcda (48) 
as measure in the wave function space. Using this measure the operators p a , p b , p c and p a are 

. d , d . d . d 



The WDW equation is 

Ft 2 Ft 2 Pi 2 _ & 1 

V(a,b,c,a) = 0. (50) 



da 2 db 2 dc 2 da 2 



The WDW equation can be completely solved by the technique of separation of variables. The 
general (bounded) solution is the superposition of wave functions 

*(a, 6, c, a) = J dk b dk c dk a A(k b , k c , k a )ip(k b , k c , k a ; a, b, c, a) , 

^(k b , k c , k a - a, b, c, a) = Ne^e^e^KUKe-") , v = sfk 2 - k 2 , k t >0, (51) 

where K iv is the modified Bessel function of imaginary index iv [24]. By properly choosing the 
normalization factor N, and fixing the gauge using the b degree of freedom, the eigenstates of 
the physical Hamiltonian with energy E = k\ read 



K***, = J^^e^e^KUKe-n . (52) 
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Note that v = k on shell. The wave functions (52) form an orthonormal basis in the Hilbert 
space with (gauge- fixed) measure d[zu] = dadcda 

(V>i, ip 2 ) = J dadcda ip*ip 2 = 5^- u 2 )5(k cl - k c2 )5{k al - k a2 ) ■ (53) 

This completes the quantization in the hybrid representation. 

Using the generating function (35) the equivalence between the two sets of (gauge off-shell) 
solutions of the WDW equation (46) and (51) can be proven. The relation between the wave 
functions in the two representations is 

#(a, b, c,a) = J d[w(u, x, y; i)] e ^(«.^^«^*%( U) Xj y - 1 ) , (54) 

where d[w(u, x, y; t)] is the measure in the space of wave functions. From Eq. (43) we know 
that d[zu(u,x,y;t)] must be 

d[w(u, x, y\ t)] = dudxdydt t~ p ~ 2 . (55) 

The unknown parameter p in Eq. (55) represents possible factor ordering ambiguities in the 
relation between the (gauge off-shell) wave functions in the two representations. This ambigu- 
ity can be fixed by requiring the equivalence of the (gauge on-shell) wave functions in the two 
representations, thus proving the physical equivalence of the representations. By substituting 
Eq. (46) and Eq. (35) in Eq. (54) and choosing p — —1/2 in Eq. (55), after some algebra 
(details are given in Appendix B) one obtains Eq. (51). 

Let us briefly discuss the correspondence between the hybrid wave functions and the clas- 
sical solutions. The oscillating regions of the wave functions correspond to the classically 
allowed regions of the configuration space. Along c and a directions the wave functions (52) 
are described by plane waves. Along the a direction the wave functions are oscillating in the 
region 

0<e" a <^. (56) 

This corresponds to the classically allowed region for the hybrid variable a. (We have chosen 
k a > for simplicity.) Indeed, from the solutions of the equations of motion we have 

< e~ a = — [cosher)]- 1 < — . (57) 

Pa ' Va 

Finally, the wave functions go like e ±iau for large values of a. With aid of Eqs. (31), Eqs. 
(32), and Eqs. (34), we find the relation between the quantum numbers ki and the classical 
parameters that characterize the behaviour of the classical solution 

k b = -2V3£ , k c = -^= , v = k. (58) 

2y 3 
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Using the previous relations it is straightforward to identify the hybrid wave functions that 
correspond to the classical cases Finally, starting from Eq. (52) quantum M-theory cos- 

mology can be investigated with aid of usual elementary quantum mechanics techniques. For 
instance, wave packets can be constructed along the lines of Refs. [25] and specific boundary 
conditions for the wave function of the universe can be selected. 

It is curious to note that the hybrid wave functions (52) are formally equivalent to wormhole 
wave functions. [See e.g. [25] and references therein.] It would be interesting to investigate 
whether this analogy is purely accidental or conceals a deeper relation between M-theory 
cosmology and wormhole physics. 



4 Invariant manifold II 



In this case the phase space is six-dimensional and the canonical action is 

S7 = J dt ap a + (fiptp + Ppp - H , 

where the Hamiltonian is 

The canonical equations of motion are 



a 



Pa 
6 ' 



r 2 ' 



P 



P± 

12 



2 2 
The previous equations are supplemented by the constraint 

2p' - Qpl + Pl + 12Q V a -^ 6/3 = 



(59) 



(60) 



(61) 



(62) 
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(63) 



Let us first discuss the case Q ^ 0. The off-shell solution of the canonical equations is 

a = a - - In [cosh («(r - r ))] - £(r - r ) , 

Pa = 2" taIlh ~ ~ 6 ^ ' 

(j) = 0o - ^ In [cosh (k(t - t ))] - x(r - r ) , 
= -tanh[«(r-ro)] + 2x, 

/3 = /3 + - In [cosh (k(t - r ))] - p(r - r ) , 

Py3 = 3/ttanh [k(t — r )] — 12p , 
where the constants of motion are related by 

k 2 + 6£ 2 - 2x 2 + 12p 2 = 2H , k ^ , (64) 
and (we choose k > for simplicity) 

3a o -0o-6/3 o = 21n(j^) , 3£ - 6p - X = ■ (65) 

On the (physical) shell x — is a degenerate trivial configuration of the system because it 
implies k = 0, £ = 0, and p = 0. 

Analogously to the invariant manifold I, we can find a maximal set of gauge invariant 
observables. In this case we expect four physical gauge invariant observables because the 
system possesses three canonical degrees of freedom. A possible choice is 

U = a - _ Pa ± ^ arC cosh ( " e -<3"-<tf-*)/^ 
Qk \\Q\ ) 

(66) 

X = I(3« + 2(3 - 0) - 3Pa + PP + 3p * arccosh (^-^-V-tvA > 

4 z4/t \ I ^ I / 
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where 



k = 



16 



{-Pa + Pf3 ~P<t>) + Q e 



2 3a-6(3-(j> 



1/2 



(67) 



The quantities (66) have been chosen such that they satisfy the canonical Poisson brackets 



i . 



The conjugate of the Hamiltonian that completes the Shanmugadhasan chart is 



T 



1 



arccosh 



K 



K 

W\' 



-(3a-6/3-0)/2 



(68) 



(69) 



Using Eqs. (60), (66) and (69) we find the relation between the gauge invariant observables 
and the constants of motion in Eqs. (63) 



a = i(3X-f/) + iln-^-, /3 = ^-Iln-^, O = l(X - U) + \ In , 



\Q\ 



\Q\ 



\Q\ 



Z = --(y + 3W/4), 



and 



W 

P = ~24> 



1 



(70) 



X = -^(V + W/4), 



K 



— {W 2 - W 2 ) + 2H 
1 2 



1/2 



(71) 



The relation between the invariant submanifold I and the invariant submanifold II can 
be made manifest in the hybrid canonical chart. For the invariant submanifold II the hybrid 
variables (a, b, c) are 

V3 



a = — — (3a — 6/3 — 0) , b = y/3(a — <p) , 



c=—(3a + 2p-<f>) , 



Pa = | {-Pa + P(i 



Pb 



(Pa + 3fty) , p c 



1 



(72) 



(tya + Pf3 + 3fty) . 



2V3 ^ ™ ' ye 4v/3 
Indeed, using the hybrid variables the constraint (62) reads (we have divided by a factor 12) 

p 2 a-pl+pl + Q 2 e- 2a = 0. (73) 

Moreover, the gauge invariant canonical variables (66) are related to the hybrid variables (72) 
by the canonical transformation 



U 



X 



1 

i 



b + ^arccosh 



K 



K 

W\ 



2^/3 



c — — arccosh ( — — e Q 



\Q\ 



T = — arccosh 

K 



K 

w 



v = Vs Pb , 

W = 2y/3 Pc , 

H = l(pl-p 2 b +p 2 c + Q 2 e- 2a ), 



(74) 
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Figure 4: Parameter space for invariant submanifold II. The physical points are represented by the 



two branches of the hyperbola p = ± ^/4£ 2 + 5k 2 /12J /5.The colored regions determine on the 
two branches of the hyperbola the different physical cases described in the text. 



where k = \Jp1 + Q 2 e 2a . Note that this transformation coincides with the section p a = Q of 
the transformation (34). 

The qualitative discussion of the physical properties of the solution proceeds similarly to 
that of the first invariant manifold. 4 The dynamics of the model is determined by k and £. 
At fixed k, the physical points of the model are determined in the plane (£, p) by the two 
branches of the hyperbola (64) 



The upper branch is characterized by an ever increasing four-dimensional effective coupling 
g, which evolves from zero at r — > — oo to infinity at r — > oo. Conversely, the lower branch is 
characterized by an ever decreasing four- dimensional effective coupling g, which starts with an 
infinite value at r — > — oo and vanishes at r — > oo. For each of the two branches we distinguish 
three different dynamical behaviours of the external geometry according to the value of £ (see 
Fig. 4). For the upper branch we have: 

4 See also Ref. [16] where eleven-dimensional cosmological solutions with and without four-form charges are 
investigated. Note that Kaloper, Kogan and Olive work in the eleven-dimensional supergravity frame whereas 
we work in the string frame. The relation between the quantities used in our paper and those used in Ref. 
[16] (denoted with a twiddle) area = a + j/2, (3 = /3 + j/2, $i = 37, and N = he^l 2 . 
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i) Ever expanding external spacetime. This happens for £ < —k/A. The ten-dimensional 
coupling constant is zero at r — > — oo and becomes infinite at r — > oo. The Hubble 
constant is always positive and vanishes at r = — oo where the external scale factor is 
also zero. At large times (r — > oo) both the Hubble constant and the external scale 
factor become infinite except for the limiting value £ = —k/A for which they reach zero 
and a finite nonzero value, respectively. The radius of the internal space is infinite at 
large times and zero, finite, and infinite for £ < -k(3/4 + 1/V3), £ = -k(3/4 + l/\/3), 
and — k(3/4 + 1/a/3) < £ < — k/4 at r — > — oo, respectively. 

m,) Ever contracting external spacetime. This happens for £ > «/4. The ten-dimensional 
coupling constant is infinite at r — > — oo and vanishes at r — > oo. The Hubble constant 
is always negative, in particular is zero at r = — oo and infinite at r = oo where the 
external scale factor is zero. At small times (r — > — oo) the external scale factor is 
infinite except for the limiting value £ = k/A for which is finite nonzero. The radius of 
the internal space is infinite at r = — oo and zero at r = oo. 

Hi) Bouncing external spacetime {—k/A < £ < k/A). In this case the radius of the external 
spacetime is zero at r — > — oo, reaches a maximum, and then decreases to zero at r — > oo. 
The Hubble constant evolves from zero to — oo passing through a positive maximum. 
The ten-dimensional coupling constant and the scale factor of the internal space evolve 
as follows 

a) For — k/A < £ < —k/12 g w evolves from zero [from a finite nonzero value for the 
limiting value £ = —k/12] to infinity and Rt& bounces from infinity to infinity; 

b) For —k/12 < £ < k/12 both g w and Rt& bounce from infinity to infinity [gio to a 
finite nonzero value for the limiting value £ = k/12] (see Fig. 5); 

c) For k/12 < £ < k(3/4 — l/Vs) gw evolves from infinity to zero and R T & bounces 
from infinity to infinity; 

d) For k(3/4 — I/a/3) < £ < «/4 both g w and i? r e evolve from infinity to zero [R T n 
to a finite nonzero value for the limiting value £ = re(3/4 — l/\/3)]. 

For the lower branch we have: 

i) Ever expanding external spacetime for £ < —k/A. The ten-dimensional coupling con- 
stant is zero at r — > — oo and becomes infinite at r — > oo. The Hubble constant is 
always positive and vanishes at r = oo where the external scale factor is infinite except 
for the limiting value £ = —k/A for which it reaches a finite nonzero value. At small 
times (r — > — oo) the Hubble constant is infinite and the external scale factor is zero. 
The radius of the internal space evolves from zero to infinity (see Fig. 6). 

ii) Ever contracting external spacetime for £ > k/A. The ten-dimensional coupling constant 
is infinite at r — > — oo and vanishes at r — > oo. Both the external scale factor and the 
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Figure 5: Graphical representation of the dynamics of the invariant submanifold II for — rc/12 < £ < 
k/12 (upper branch) as function of the gauge proper time t g and of the cosmic proper time t c . The 
external scale factor is vanishing at small times, reaches a maximum, and then decreases to zero at 
large times. Both the internal scale factor and the coupling constant g\o are infinite at small times, 
decrease to a minimum value, and then increase again to an infinite value for large times. (See Fig. 
2 for the key to the plotted lines.) 




Figure 6: Graphical representation of the dynamics of the invariant submanifold II for £ < — k/4 
(lower branch) as function of the gauge proper time t g and of the cosmic proper time t c . The 
scale factor of the external spacetime, the scale factor of the internal space and the ten-dimensional 
coupling constant gio are ever increasing. The coupling constant g is ever decreasing. (See Fig. 2 for 
the key to the plotted lines.) 
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Hubble constant (always negative) vanish at large times. At small times the radius 
of the external spacetime is infinite and H(—oo) = — oo except for the limiting value 
£ = re/4 for which the radius of the external spacetime and H are finite nonzero and 
zero, respectively. In this case H evolves through a negative minimum before going to 
zero at large times. The radius of the internal space is infinite for r = — oo and infinite, 
finite nonzero, and zero at r = oo for re/4 < £ < re(3/4 + I/a/3), £ = re(3/4 + l/y/3), 
and £ > re(3/4 + 1 / a/3), respectively. 

Hi) Bouncing external spacetime for —re/4 < £ < re/4. In this case the radius of the external 
spacetime is zero at r — > — oo, reaches a maximum, and then decreases to a zero value at 
t — > oo. The Hubble constant evolves from infinite to zero through a negative minimum. 
The ten-dimensional coupling constant and the scale factor of the internal space evolve 
as follows 

a) For —re/4 < £ < — re(3/4 — l/\/3) <?io evolves from zero to infinity and R T & evolves 
from zero [from a finite nonzero value for the limiting value £ = — re(3/4 — l/\/3)] 
to infinity; 

b) For — re(3/4 — l/\/3) < £ < —re/12 gio evolves from zero [from a finite nonzero 
value for the limiting value £ = —re/12] to infinity and R T <> bounces from infinity 
to infinity; 

c) For —re/12 < £ < re/12 both g w and Rt& bounce from infinity to infinity [qiq to a 
finite nonzero value for £ = re/ 12]; 

d) For re/12 < £ < re/4 gio evolves from infinity to zero and Rtb bounces from infinity 
to infinity. 

From the physical point of view cases i) (ever expanding external spacetime) are not viable 
because both the eleven dimension and the internal space always decompactify at large times. 
Bouncing external spacetimes can instead give a more realistic description of our universe. 
In particular, case iii)-d) of the upper branch is characterized by an ever decreasing ten- 
dimensional coupling and an ever decreasing internal scale factor. This property guarantees 
that a large-size external spacetime is weak-coupled and its internal space is compactified. 
Moreover, a bouncing external spacetime is always characterized by early accelerated expan- 
sion. Models with bouncing g w and R T a are probably also viable from the physical point of 
view, albeit they might require some kind of fine tuning of the free parameters to actually 
describe a realistic universe. 

4.1 Quantization 

The quantization of the system proceeds analogously to the case I. (See also the footnote in 
Subsection 3.2.) Let us start again with the Shanmugadhasan representation. By imposing 
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the gauge fixing (39) we obtain the effective action 

S cS = Jdt [UV + XW- H eS ] , (75) 

with vanishing effective Hamiltonian H c r = —H. We work in the representation in which 
V, W are differential operators with eigenvalues v and w, respectively. The Hilbert space 
measure is 

d[w] = dudx , (76) 

and V, W are 

v = > w = > ( 77 ) 
ou ox 

An orthonormal basis in the Hilbert space is given by the set of eigenfunctions of the observ- 
ables (77), namely 

V(u,x,y) = -l-e i ( wl+wx ) (78) 

In the hybrid representation the quantization is formally identical to that of case I. So we 
give here only the result. The Hamiltonian eigenstates with energy E — k\ are 



i>(k b , k c ; a, c) = ^^" e^KUQe-") , v=y/tf-k* c , k t > , (79) 

The set of wave functions (79) forms an orthonormal basis w.r.t. gauge- fixed Hilbert space 
with measure d[zu] = dadc. Starting from Eqs. (79) the quantum mechanics of the invariant 
submanifold II can be investigated. 

Finally, let us discuss the degenerate case S Inn (Q — 0). The general (Kasner) solution is 
(see also [16] and footnote 4) 

a = a + y(r - r ) , = O - y (r - r ) , (3 = (3 + ^(r - r ) , (80) 

where p a , p^, and pp are constants related by the condition 

2p 2 a -6pl+p} = 2H. (81) 

Note that Eqs. (63) reduce to Eqs. (80) with k — and a redefinition of the integration 
constants. The on-shell dynamics is identical to the dynamics of a Klein-Gordon particle 
moving in a three-dimensional Minkowski space with timelike coordinate tp/y/E and spacelike 
coordinates a / \/2 and (3. In this case the gauge and hybrid canonical variables coincide and 
the quantization of the system is straightforward. Classically, the SVn// model describes a 
spacetime with expanding, contracting, and constant internal (external) space for pp > 
{p a > 0), p/3 < (p a < 0), and pp = {p a = 0), respectively. The sign of p<p determines the 
strong and weak coupling regions of the model. In particular, a geometry with contracting 
internal space (pp < 0) and expanding external spacetime evolves from a strong region to a 
weak region for p^ > p a + pp > and —pp/2 < p a < —5pp/2 (see Fig. 7). 
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Figure 7: Parameter space for the degenerate case Smii- The portion of the hyperbola in the 
colored region represents geometries with contracting internal space (pp < 0) and expanding external 
spacetime which evolve from a strong coupling regime (g — > oo, 510 — > 00) at r = — 00 to a weak 
coupling regime (g « 0, 510 ~ 0) at r = 00. 

5 Conclusions 

In this paper we have analysed a spatially flat, four-dimensional cosmo logical model derived 
from the M-theory effective action [Eq. (1)]. The eleven-dimensional metric is first compact- 
ified on a one-dimensional circle to obtain the type IIA superstring effective action and then 
on a six-torus to obtain the effective four-dimensional theory. In our investigation we con- 
centrated the attention on the boundary of the physical phase space of the theory which is 
constituted by two invariant submanifolds, where either the axion field (from the NS/NS sec- 
tor) is negligible or the four-form R/R field strength is irrelevant. In our discussion we have 
heavily employed the canonical formalism. This approach makes the analysis of the features 
of the classical solution extremely simple and allows a straightforward quantization of the 
theory. 

We found regions in the moduli space where a four-dimensional FRW universe evolves 
from a strong coupling regime towards a weak coupling regime, both internal six-volume and 
eleven-dimension contracting. In some cases, the dynamics is characterized by an early accel- 
erated (inflationary) expansion with the spacetime eventually approaching a standard FRW 
decelerated expansion. Such scenario happens when the pseudo-scalar axion field dominates 
and the R/R field strength is negligible. (See also Ref. [15].) When the axion is subdominant 
w.r.t. the R/R field, different cosmological scenarios appear on the scene. In this case an ever 
expanding external spacetime is characterized at large times by an infinite ten-dimensional 
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coupling constant (the size of eleven dimension). This seems to rule out this model from a 
physical point of view. On the contrary, for bouncing external spacetimes, where the external 
dimensions recollapse after initially expanding, geometries with early accelerated expansion, 
decreasing ten-dimensional coupling constant and decreasing internal space are found. 

The quantization of the two invariant submanifolds can be performed exactly and the 
Hilbert space of states is derived. In the quantum framework our analysis allows to iden- 
tify the quantum states that correspond to the different classical behaviours. In the hybrid 
representation we have identified regions in the space of parameters where the wave function 
of the universe is either oscillating or exponentially decaying. These regions are determined 
by the inverse exponential function of the four-dimensional (unshifted) dilaton, i.e., by the 
four-dimensional string coupling, and correspond to classically allowed and classically forbid- 
den regions, respectively. Starting from the Hilbert space of states, the quantum mechanics 
of M-theory cosmology can be constructed with aid of usual elementary quantum mechanics 
techniques. 

The results summarized in the preceeding paragraphs strengthen and extend previous 
results obtained within different approaches [13]-[17] . Although the models investigated in 
this paper represent a severe, yet consistent, truncation of the original eleven-dimensional 
action Eq. (1), new interesting physical output can be extracted. In particular, the quantum 
mechanics of M-theory that has been presented here may constitute the basis for further 
insights in high-energy (quantum) M-theory cosmology. If M-theory is the ultimate theory of 
quantum gravity, M-theory quantum cosmology is a subject which is certainly worth exploring. 
From this point of view related topics of investigation which deserve attention are the study 
of inhomogenous perturbations of the fields and a complete analysis of FRW models in the 
full phase space of the theory, when all form fields are excited. In the latter D amour 

and Henneaux have recently observed [12], the dynamics may be qualitatively different and 
lead to potentially new theoretical and observational effects. 
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Appendix A. Comparison with notations of BCLN. 

In this Appendix we compare our notations with those of Ref. [17]. This allows a straightfor- 
ward interpretation of our results in the formalism of BCLN. 
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Let us first consider the invariant submanifold I. Using Eqs. (16) and (23), and Eq. (26) 
of Ref. [17] the BCLN variables 

x bcln, Ubcln an d Qbcln read 



Pa 

%BCLN ' 



respectively. Moreover we have the important relation 



7/ - PP O - 2 (— \ p - 2 ( 3q +^) (R2) 

VQp<f> \P<t>J 



4 



^ BCLN — 1 X BCLN V BCLN ^BCLN + 2^ ~ ^ ' (^3) 

P<f> 



Differentiating Eqs. (82) w.r.t. (j> an d recalling Eq. (23) we have (on shell) 



Qbcln{%bcln "I" V^3) , 



— QbCLnD BCLN 1 (^4) 



' BCLN 

dy 

BCLN 



d>T BCLN 

dQ BCLN _ / 2 , 2 , /o \ 

— z,±L BCLN yx BCLN T SJbcln ' V OX BCLN) 1 

(XT, 



BCLN 



where we have used Eq. (83) and the equations of motion (16). Eqs. (84) coincide with Eqs. 
(28)-(29) of Ref. [17] for the invariant submanifold I. Since Q BC ln 

> from Eq. (83) the 

physical region in the phase space is 

x 1cln + Ubcln ^ 1 > (85) 

thus recovering the BCLN result. Finally, from (84) we obtain 
d 



dr E 



I %BCLN ~\~ BCLN 

+ y/3) = Vt BCLN (x BCLN + V2y 

BCLN 

+ V3 . (86) 

' BCLN 

Equation (86) coincides with Eq. (33) of Ref. [17] for z = 0. 

BLCN present the analitical solution for the invariant submanifold I (see Eq. (35) of Ref. 
[17]). Using Eq. (82) and the solutions of the Eqs. of motion (24) it is straightforward to 
obtain Eq. (35) of Ref. [17] where 

vo _ pp (g7) 



Let us now consider the invariant submanifold II. Using Eqs. (61) and (23), and Eq. (26) 
of Ref. [17] the BCLN variables x bcln, y bcln and z BC ln read 

r - Pa y z - o ( 9.Y e^-i-V 

^ bcln — pr i Ubcln — /— ; ^bcln — ^ I I c i 

V3ps VQps \P<t>J 



27 



respectively. Using Eqs. (88) we find 



4 



&>BCLN — 1 X BCLN V BCLN Z BCLN + 2 ^ ~ ^ ' (^9) 

P<f> 

thus Vt BCLN is vanishing on-shell as expected for the invariant submanifold II. Since z BCLN > 
the physical region in the phase space is 

X BCLN + VbCLN — 1 ' (^O) 

where the equality holds iff Q = 0. 

Differentiating the previous equations w.r.t. and recalling Eq. (23) we have 

dXgCLN Z BCLN 
dT BC LN 2 

dy 

BCLN Z BCLN 



dT BC LN 2 
BCLN 



1 BCLN 

dz 



X BCLN v3J , 

(j/bcln + V0) , (91) 



Z BCLN {^BCLN 1 V^6|/sCL]V + V^EbCLJV") 



Eqs. (91) coincide with Eqs. (28)-(29) of Ref. [17] for the invariant submanifold II. [Remember 
Eq. (89)]. Finally we have 

+ V2y BCLN + V3 = —L— (3p^ - p a - pp) , (92) 
V3p<f> 



x 



BCLN 



and by differentiating Eq. (92) 

~l [XgCLN "i" 

V2y 

BCLN + VS) = (x BCLN + V2y 

BCLN ~ 73 . (93) 
dr BCLN v ' 2 v y 

Equation (93) coincides with Eq. (33) of Ref. [17] when Vt BCLN = 0. 

BLCN present also the analitical solution for the second invariant submanifold (Eq. (34) 
of Ref. [17]). Again, it is straightforward to obtain the BCLN result where 

^±^ = V2^±P. (94) 
This completes the comparison to the BCLN variables. 



Appendix B. Proof of the equivalence between gauge in- 
variant and hybrid quantum representations. 

Let us consider Eq. (54). By substituting Eq. (46) and Eq. (35) the wave function in the 
(a, b, c, a) representation can be written 

$(a,b,c,a) = N J dtr^hiu.t-^hix.t-cjhiy.t-a.a), (95) 
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where 



Ii(u,t;b) = J duexp 
J 2 (a;,t;c) = J dxexp 

I 3 (y,t;a,a) = J dyexp 



tuv + i 



2t 



ixw — i- 



. (2VSx - c) 



2t 



xyz — i 



arccosh 2 ^! + e 2a (y — a)' 



2t 



(96) 
(97) 

(98) 
(99) 



and iV is a normalization factor. The evaluation of the integrals I\ and J 2 is immediate. The 
result is 5 



1 1 = \ft exp 

1 2 = Vi exp 



/ bv tv 2 \ 

( cw tw 2 
\2V3 + 



(100) 
(101) 



Substituting the previous results in Eq. (95) and using the new variable y = y — a, Eq. (95) 
is cast in the form 



\I>(a, 6, c, a) = N exp 
where 



bv cw 

— 1= H y= + <JZ 

2V3 , 



dye iVz J 



dtt'^exp 



w 2 — 4v 2 f- arccosh y/1 + e 2a y 2 

S = — » C = 2V3«- 



, (102) 



(103) 



12 ■ " ' Vw 2 - 4v 2 

From Ref. [24] (Eq. 8.421 No. 8, p. 966) the integral in dt can be evaluated. Using the new 
integration variable x — arccosh y/1 + e 2a y 2 Eq. (102) is cast in the form 

bv cw 



fy(a,b,c,a) = N exp 



V3 2V3 , 



J dx cosh XX P H { p 



\Jw 2 — 4t> 2 " 



2^ 



(104) 
(105) 



Choosing p = -1/2 and recalling that H W 1/2 (z) = ^/2/irze iz (see Ref. [24] Eq. 8.469, No. 6, 
p. 978) we have 



\&(a, 6, c, a) = iVexp 



bv cw x 
—j= H j= + oz 



± 



duj 



±1 

arcsinh uj+ize a u> 



(106) 



3 Here and throughout the Appendix neglect overall constant factors that are reabsorbed in N. 
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where u = sinhx. Finally, from Eq. 3.483, p. 387, of Ref. [24] we find 

' bv cw 



ty(a,b,c,a) = N exp 
= iVexp 

where 



. I ,- H 7= + crz 

\V3 2v/3 , 

. / bv cw 

Vv3 2V3 , 



(ze- a K v+l (ze~ a ) + ze^K^ze" 1 )) 
K(ze- a ), (107) 



v 



\Jw 2 — 4v 2 
2^/3 



(108) 



and we have used Eq. 8.471, p. 979, of Ref. [24]. Setting v = ±y/3k b , w = ±2y^3k c , and 
z = ±k a (kb, k c , k a > 0) we obtain Eq. (51). 
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